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Abstract 


The focus of operations research is the existence of a problem that requires making an appropriate decision that helps 
reduce risk and achieves a good level of performance. Operations research methods depend on formulating realistic 
issues through mathematical models consisting of a goal function and constraints, and the optimal solution is the ideal 
decision, despite the multiplicity of these methods. However, we encounter many complex issues that cannot be 
represented mathematically, or many issues that cannot be studied directly. Here comes the importance of the 
simulation process in all branches of science, as it depends on applying the study to systems similar to real systems 
and then projecting this. The results if they fit on the real system. So simulation is the process of building, testing, and 
running models that simulate complex phenomena or systems using specific mathematical models. The simulation 
process depends on generating a series of random numbers subject to a regular probability distribution in the field [0, 
1], and then converting these random numbers into random variables subject to the distribution law. Probability, 
according to which the system to be simulated operates, using appropriate techniques for both the probability density 
function and the cumulative distribution function. Classical studies have provided many techniques that are used 
during the simulation process, and to keep pace with the great scientific development witnessed by our contemporary 
world, we found that a new vision must be presented for this. Techniques A vision based on the concepts of 
neutrosophics, the science founded by the American mathematical philosopher Florentin Smarandache. The year 1995, 
in which new concepts of probabilities and probability distributions are used, as we presented in previous research 
some techniques from a neutrosophic perspective, and as an extension of what we presented previously, we present in 
this research a neutrosophic vision of the Box and Muller technique used to generate random variables that follow a 
normal distribution. 


Keywords: Simulation; Neutrosophic simulation; Neutrosophic normal distribution; Generation of neutrosophic 
random variables that follow the normal distribution; Box and Muller technique 


1. Introduction 
In the world of rapid changes, there has emerged an urgent need to make a rational decision based on quantitative 
methods and methods that limit the risk ratio, especially if the decisions are crucial and the decision issues are huge 
and complex, which was the reason for the emergence of the science of operations research, the science that is the 
applied side of mathematics and helped improve the performance of many systems that used its methods in their 
workflow. One of the most important methods of operations research is the simulation method, which is a numerical 
technique used to carry out tests on a numerical computer and includes mathematical relationships that interact with 
each other to describe the behavior and structure of a system in the real world over a period of time without any risk 
or procedure. Experiments on the real world and allow the study of complex issues that are difficult to solve using 
other operations research methods. We know that the basis of the simulation process is generating a series of random 
numbers subject to a regular probability distribution in the range [0,1], and then converting these random numbers 
into random variables that follow the distribution. The probability according to which the system to be simulated 
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works. In classical logic, there are many techniques that can be used in the conversion process, which result in random 
variables that follow irregular probability distributions, to keep pace with the great scientific development that our 
contemporary world is witnessing and after the emergence of neutrosophic logic and the accurate results presented by 
research and studies. presented by researchers and those interested in this logic in most fields of science, including [1- 
13], and given the great importance of the simulation method, we presented in previous research [14-18] a 
neutrosophical vision of some of the techniques used to generate random variables that follow important probability 
distributions, such as the inverse transformation technique. Which can be used for all probability distributions whose 
cumulative distribution function has an inverse function, such as the uniform distribution over the domain [a, b], the 
exponential distribution, the Weibull distribution, and others. We also presented a neutrosophical vision of the 
rejection and acceptance technique and how to use it to generate random variables that follow the beta distribution. 
We present in this research is a neutrosophical study of the Box and Muller technique used to generate random 
variables that follow complex probability distributions. The aforementioned techniques cannot be used in the process 
of generating random variables that follow a normal distribution. 


2. Discussion: 

2.1. Box and Muller technique for generating neutrosophic random variables that follow a normal 
distribution: 

Within this technique, we use many mathematical relationships in order to be able to represent some complex 
probability density functions in a more simplified way so that random variables can be generated using known 
techniques. We know that the simulation process depends on generating a series of random numbers that follow a 
uniform distribution over the field [0,1], and using the appropriate transformation method we convert these random 
numbers into random variables that follow the probability distribution according to which the system to be simulated 
operates, based on which to generate random variables. Neutrosophic, we generate a series of random numbers that 
follow a uniform distribution over the field [0,1]. From this series we can obtain neutrosophic random variables that 
follow the probability distribution according to which the system to be simulated operates. Here we distinguish three 
cases: The first case: Neutrosophic random numbers and probability distribution in the classical form The second case: 
classical random numbers and the probability distribution is in the neutrosophic form. The third case: neutrosophic 
random numbers and probability distribution in neutrosophic form. The normal distribution is one of the most widely 
used probability distributions, but from the definition of its probability density function, we find that it is not possible 
to use the neutrosophic study presented in research [15,17] for the inverse transformation method and the rejection 
and acceptance method. Rather, we need a special technique for probability distributions in which the probability 
density function is a complex function, such as the Box-Muller technique. 

2.2. neutrosophical study of the Box and Muller technique for generating random numbers following a 

general normal distribution: 

2.2.1. | This technique was presented according to classical logic in many references, including references [19,20], 
and its scientific basis was to derive two simple probability distributions from the general normal distribution 
after converting it to the standard normal distribution. 

Dedede Neutrosophic study: 

a. In Reference [21], the neutrosophic probability density function for a random variable following a normal 
distribution is defined by the following relation: 


: 1 (2) 
Xv~N, ,0, = ex: 20N 1 
nv~Ny (tw, On“) Sane (1) 
b. We convert it to the standard normal distribution using the following relation: 
a Sed Ly) 
iy =—— (2) 
On 
c. Then the neutrosophic probability density function for the standard normal distribution is written as follows: 
1 _4N 
Zy) =——=e 2 3 
f (Zn) pr (3) 


d. Using the relationship (3), we form two independent neutrosophic variables Z,, and Zz, subject to the 
standard normal distribution as follows: 
fZw)=aee F fy) =e 
= (2 —I———| eZ 
1N - 2N = 


Forming their joint distribution, we find: 
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1 lz? +22y) 
f Zin Zon) = f Zin)- f Zon) = ee Cee ie 
e. We move from Cartesian coordinates to polar coordinates using transformations: 
Zin =TyCOSO = Zoy = Tysind 
Substituting this transformation and its extensions |J| = rj), we obtain the following joint distribution in terms of the 
polar variables ry, 0: 


1 2 
fw, 9) = ae Bu i 


From it we obtain the following two single neutrosophic distributions for Ty, 0: 


sd iy pmaes oe se sae al 
ro =| Pein Odry == | Ty e BON ae —e2 , Se 


f. The first distribution is a uniform distribution over the domain [0,2p] defined by the following relation: 
il 
é)=— 4 
f@=— 
In the same way we find: 
21 


TN 
Tye 2 (27 TH 
fw =] f@y,@)de = I d@ =tTye 2 
0 2m Jo 


It is a Rayleigh distribution defined on the field [0, +00[. To convert this distribution to an exponential distribution, 
we perform the following transformation: 


: dry 1 
Uy = Ty = duy = SINGIN oe Fae Ore 
Substituting the Rayleigh distribution, we get: 
_un 1 1 _un 
fun) =™me 2. Tee 4 x 
g. The second distribution is the exponential distribution whose parameter is J and is defined by the following 


relation: 


1 uy 
fly)=5e 2 () 
2.2.3. From the above, the Box and Miiller technique is formulated in the following steps: 
1) We generate two random numbers R, and R, that follow a uniform distribution over the range [0,1]. 
2) We convert the two random numbers into two neutrosophic random numbers, R,y and Ry, see [14]. 
3) We use the inverse transformation technique to transform the first number R,,y into a random variable that 


follows an exponential distribution whose parameter A = ~ using the following relation: 


—InRiy 
Ea ee 
We get: 
al 
rh = —2InR,y = Ty = (—2InR,y)2 


4) Weuse the inverse transformation technique to transform the second number Ry into a random variable that 
follows a uniform distribution over the field [0,277] using the following relation: 
Xiy = (b—a)Riy ta 
We get: 
Oy = 27R2y 
We substitute the following two relations: x,y = Tycos@ ,xX2y = TysinO we get: 


1 
Xin = Care cos(2mRzy) (6) 


Xan = (—2InR,y)2sin(2mRzy) (7) 
2.2.4. The Box and Muller technique for generating random variables follows the general normal 
distribution N (Hy, a), defined by the following relation: 


Pe ee ae 
rae 


We transform random numbers subject to the standard normal distribution by the relation: 
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Xy = Oy-Zy + My 
Thus, we get the following generation relation: 


1 
Xy = Oy.[—InRzn]2 cos(27Ryy) + My (8) 
Practical example: 
Using the Box and Muller technique, find two random variables that follow a normal distribution with uw = [15,17] 
and o = [2,3]. 
Using the mean square method, we generate two regular random numbers in the range [0,1]. 
The mean square method is given by the following relation: 
Riz, = Mid[R?] ;i = 0,1,2,3, ... 
Where Mid represents the middle four levels of R?, and R; is chosen, any random fractional number consisting of four 
places (called a seed) that does not contain a zero in any of its four places, we get: 
R, = 0.1234 ,R, = 0.5227 
We convert these random numbers into neutrosophic random numbers and here we distinguish three states of the field 
[0,1] with a margin of indeterminacy, in the three states we have € € [0,n] andO <n<1 
First case: [0 + €, 1] indeterminacy at the minimum of the range. We substitute in the relation: 


R,-€ 
Riy = = 
Second case: [0,1 + €] indeterminacy at the upper limit of the range. We substitute in the relation: 
Pin = pace 
IN 1+¢ 
Third case: [0 + €,1 + e€] Indeterminacy in the upper and lower limits of the field. We substitute in the relation: 
Riy = R, — (s 


Let us take € € [0,0.02] then we get: 
In the first case: 


R,—e€ _ 0.1234—[0,0.02] _ [0.1034,0.1234] 


Ry = —— = — = = 0103 4,0.1259 
WW 1-8 1 — [0,0.02] [1,0.98] 

—_ R,-€ — 0.5227 —[0,0.02] _ [0.5027,0.5227] _ 10.5027,05334] 
aN 1-e ° .©}©1-—[0,0.02] — [1,0.98] — a 


We replace the relationship (8) 


1 
Xiy = Oy-[—INRzy |? cos(27Ryy) + Ly 
A: 
x1y = [2,3][—In[0.5027,0.5334]]?. cos(27)[0.1034,0.1259] + [15,17] 
Xy = [16.6584,19.3782] 
We compensate with the relationship (9) 


4: 
Xeon = On: [—lnR2y]2 sin(27R,y) + Ly 
al: 


Xoy = [2,3][—In[0.5027,0.5334]]?. sin(27)[0.1034,0.1259] + [15,17] 
Xoy = [15.0187,17.0328] 

Thus, we obtained two neutrosophic random variables that follow a normal distribution with 
pe = [15,17] and o = [2,3]. 
In the same way, we can obtain the random variables according to the second and third cases. 
3. Conclusion and results: 
Due to the great importance of the simulation method, and in order to obtain more accurate results that complement 
what we presented in previous research from the neutrosophic study of generating neutrosophic random numbers and 
then converting them into random variables that follow the probability distribution with which the system to be 
simulated operates, and since the reverse transformation method and the rejection and acceptance method do not they 
can be used for probability distributions defined by a complex probability density function whose cumulative 
distribution function is impossible to obtain. In this research, we presented a neutrosophic study of the Box and Muller 
technique used to treat such distributions using mathematical relationships that enable us to derive easy probability 
distributions in which the neutrosophic inverse transformation can be used to generate neutrosophic random variables. 
Then we make appropriate substitutions through which we obtain neutrosophic random variables that follow 
probability distributions that have a complex probability density function. To clarify the algorithm, we used the normal 
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distribution, which is one of the most important and most widely used probability distributions in practical life. The 
result of the study was mathematical relationships that can be used directly to generate the required random variables. 
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